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Abstract 

We derive a microscopic relativistic point-coupling model of nuclear many-body dy- 
namics constrained by in-medium QCD sum rules and chiral symmetry. The effective 
Lagrangian is characterized by density dependent coupling strengths, determined 
by chiral one- and two-pion exchange and by QCD sum rule constraints for the 
large isoscalar nucleon self-energies that arise through changes of the quark con- 
densate and the quark density at finite baryon density. This approach is tested in 
the analysis of the equations of state for symmetric and asymmetric nuclear matter, 
and of bulk and single-nucleon properties of finite nuclei. In comparison with purely 
phenomenological mean-field approaches, the built-in QCD constraints and the ex- 
plicit treatment of pion exchange restrict the freedom in adjusting parameters and 
functional forms of density dependent couplings. It is shown that chiral (two-pion 
exchange) fluctuations play a prominent role for nuclear binding and saturation, 
whereas strong scalar and vector fields of about equal magnitude and opposite sign, 
induced by changes of the QCD vacuum in the presence of baryonic matter, generate 
the large effective spin-orbit potential in finite nuclei. 
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1 Introduction 



Concepts of effective field tlieory and density functional theory are at the 
basis of many successful nuclear structure models. Quantum Hadrodynam- 
ics (QHD), in particular, presents a field theoretical framework of Lorentz- 
covariant, meson- nucleon or point-coupling models of nuclear dynamics, that 
have been successfully employed in studies of a variety of nuclear phenom- 
ena [1]. Structure models based on the relativistic mean- field (RMF) approx- 
imation of QHD have been applied to describe spherical and deformed nuclei 
all over the periodic table [2]. More recently, this relativistic framework has 
also been used in studies of the structure of exotic nuclei with extreme isospin 
and close to the particle drip lines. 

The effective Lagrangians of QHD consist of known long-range interactions 
constrained by symmetries and a set of generic short-range interactions. The 
most successful QHD models are phenomenological, with parameters adjusted 
to reproduce the nuclear matter equation of state and a set of global prop- 
erties of spherical closed-shell nuclei. On the other hand, the description of 
nuclear many-body dynamics should ultimately be constrained by the under- 
lying theory of the strong interactions - Quantum Chromodynamics (QCD). 
It is desirable to have a more explicit connection between low-energy QCD 
and nuclear phenomenology. 

As pointed out in Refs. [1,3,4], the effective hadronic Lagrangians of QHD and 
related models are consistent with the symmetries of QCD: Lorentz invariance, 
parity invariance, electromagnetic gauge invariance, isospin and chiral sym- 
metry. However, QHD calculations do not include pions explicitly. Effects of 
two-pion exchange are treated implicitly through a phenomenological scalar, 
isoscalar mean-field. Using the linear cx-model Lagrangian to construct a dy- 
namical model for the s-wave vrvr scattering amplitude, it has been shown 
that the exchange of a light scalar boson can be viewed as an approximate 
description of the exchange of correlated s-wave pion pairs [5]. While QCD 
symmetries constrain the effective QHD Lagrangians by restricting the form 
of possible interaction terms [4], the empirical data set of bulk and single- 
particle properties of finite nuclei can only determine six or seven parameters 
in the general expansion of the effective Lagrangian in powers of the fields and 
their derivatives [6] . Such a general expansion is guided by the "naive dimen- 
sional analysis" (NDA) [7,8,9,10]. NDA tests the coefficients of the expansion 
for "naturalness" , i.e. this method permits to control the orders of magnitude 
of the coupling constants. NDA can exclude some interaction terms because 
their couplings would be "unnatural" , but it cannot determine the model pa- 
rameters at the level of accuracy required for a quantitative analysis of nuclear 
structure data. 
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The success of the relativistic nuclear mean-field phenomenology has been 
attributed to the large Lorentz scalar and four-vector nucleon self-energies, 
which are at the basis of all QHD models [3]. Even though the presence of 
large scalar and vector potentials cannot be directly verified experimentally, 
there is strong evidence, in particular from nuclear matter saturation and spin- 
orbit splittings in finite nuclei, that the magnitude of each of these potentials 
is of the order of several hundred MeV in the nuclear interior. Studies based on 
finite-density QCD sum rules have shown how such large scalar and vector nu- 
cleon self-energies arise naturally from in-medium changes of the scalar quark 
condensate and the quark density [11,12,13]. In applications of QHD mod- 
els to finite nuclei, however, the scalar and vector nucleon self-energies are 
treated in a purely phenomenological way. In some models they arise from the 
exchange of "sigma" and "omega" mesons, with adjustable strength parame- 
ters determined by the nuclear matter saturation and nuclear structure data. 
Additional nonlinear terms, or density-dependent coupling constants adjusted 
to Dirac-Brueckner self-energies in nuclear matter, are necessary in order to 
obtain a satisfactory nuclear matter equation of state and to reproduce the 
empirical bulk properties of finite nuclei. 

The approach we adopt in the present work is based on the following conjec- 
tures that are believed to establish links between QCD, its symmetry breaking 
pattern, and the nuclear many-body problem: 

• The nuclear ground state is characterized by large scalar and vector fields of 
approximately equal magnitude and opposite sign. These fields have their 
origin in the in-medium changes of the scalar quark condensate (the chiral 
condensate) and of the quark density. 

• Nuclear binding and saturation arise prominently from chiral (pionic) fiuc- 
tuations (reminiscent of van der Waals forces) in combination with Pauli 
blocking effects, superimposed on the condensate background fields and cal- 
culated according to the rules of in-medium chiral perturbation theory. 

The first item has a direct connection with QCD sum rules at finite den- 
sity. The second item is motivated by the observation that, at nuclear matter 
saturation density, the nucleon Fermi momentum kf and the pion mass 
represent comparable scales. We therefore include pions as explicit degrees 
of freedom in the description of nuclear many-body dynamics. Hammer and 
Furnsthal [14] pointed out earlier that an effective field theory of dilute Fermi 
systems must include pions for a proper description of the long-range physics. 
We support this statement by explicit calculations. While it is well known that 
the lowest order Fock terms from pion exchange are relatively small, the pri- 
mary importance of chiral dynamics in nuclear systems is manifest in two-pion 
exchange processes combined with Pauli blocking effects which lead beyond 
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standard mean-field approximations. 



The broad spectral distribution of intermediate tt — tt pairs witfi = 0+, 
exclianged between nucleons, has large low- mass components starting at 2m j^. 
Hence their propagation in the presence of the nuclear Fermi sea is resolved at 
the Fermi momenta A;/ — 2m^ characteristic of nuclear matter. Guided by such 
considerations the nuclear matter problem has been investigated in Ref. [15] 
within the framework of chiral effective field theory. The calculations have been 
performed using in-medium chiral perturbation theory to three-loop order 
and incorporate the one-pion exchange Fock term, iterated one-pion exchange 
and irreducible two-pion exchange. The resulting nuclear matter equation of 
state is expressed as an expansion in powers of the Fermi momentum kf. The 
expansion coefficients are functions of kf/rriT^, the dimensionless ratio of the 
two relevant small scales. The calculation involves one single momentum space 
cutoff A which encodes NN-dynamics at short distances. This cutoff scale is the 
only free parameter. Its value can be adjusted to the energy per particle E{kfo) 
at saturation density. With A ~ 0.65 GeV adjusted to E{kfo) = —15.3 MeV, 
the calculated equation of state gives the saturation density po = 0.178 fm~^, 
the compression modulus Kq = 255 MeV, and the asymmetry energy A{kfQ) = 
33.8 MeV at saturation density. The corresponding momentum- and density- 
dependent single-particle potential U{p,kf) of nucleons in isospin-symmetric 
nuclear matter has been calculated in Ref. [16]. In particular, it has been shown 
that, at equilibrium nuclear matter density, chiral one- and two-pion exchanges 
generate an attractive nuclear mean-field for nucleons at the bottom of the 
Fermi sea, with a depth of U (0, fc/o) = —53.2 MeV, in good agreement with the 
depth of the empirical optical-model potential deduced by extrapolation from 
elastic nucleon-nucleus scattering data, and with the average nuclear potentials 
employed in shell-model and mean-field calculations for finite nuclei. 

In Ref. [17] we have outlined an approach to the nuclear many-body problem, 
both for nuclear matter and finite nuclei, which emphasizes the connection 
with the condensate structure of the QCD ground state and pion dynamics 
based on spontaneous chiral symmetry breaking. A relativistic point coupling 
model has been used, with density dependent couplings governed by scales 
of low-energy QCD. We have shown that chiral fluctuations play a prominent 
role for nuclear binding and saturation, whereas strong scalar and vector fields 
of equal magnitude and opposite sign, induced by changes of the QCD vacuum 
in the presence of baryonic matter, drive the large spin-orbit splitting in finite 
nuclei. The magnitude of these condensate background nucleon self-energies 
has been determined by the leading order in-medium QCD sum rules. 

In a different implementation, a model emphasizing the same connections 
to QCD symmetries and structure, was introduced and applied at the one- 
baryon-loop level to finite nuclei and nuclear matter in Ref. [18]. 
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In this work we expand our previous studies and present a generalized and con- 
sistent version of the relativistic point-couphng model introduced in Ref. [17]. 
For instance, in the previous work following Ref. [15], it was still assumed 
that at nuclear matter saturation density the ratio of the scalar and vector 
condensate potentials, arising through changes in the quark condensate and 
the quark density at finite baryon density, equals —1 and that the nuclear 
binding and saturation result entirely from pionic (chiral) fluctuations. This 
approximation is relaxed in the present generalized approach and the effects of 
the condensate background fields at finite density are included already at the 
nuclear matter level. Eventually we will demonstrate that with a small num- 
ber of parameters, determined by the underlying microscopic dynamics, and 
with minor fine-tuning of higher order corrections, it is possible to describe 
nuclear matter, as well as properties of finite nuclei, on a level comparable to 
phenomenological relativistic mean-field models. 

In Sec. 2 we present the formalism and develop the relativistic point-coupling 
model constrained by QCD and chiral symmetry. An effective interaction, 
characterized by density dependent coupling parameters of the local four- 
point interaction terms, is derived in Sec. 3 and adjusted to the equation 
of state for symmetric and asymmetric nuclear matter. In Sec. 4 the model 
is employed in self-consistent calculations of ground-state properties of light 
and medium-heavy N ^ Z nuclei. The calculated bulk properties and single- 
nucleon spectra are compared with experimental data and with results of the 
standard boson-exchange relativistic mean-field model. Summary and conclu- 
sions are presented in Sec. 5. 



2 Point-coupling model 

2. 1 Lagrangian 

The relativistic point-coupling Lagrangian is built from basic densities and 
currents bilinear in the Dirac spinor field ip of the nucleon: 

ilOrTi/j , Ore{l,Ti} , r G {1,7^,75,757^,(7/,^} . (1) 

Here Tj are the isospin Pauli matrices and F generically denotes the Dirac ma- 
trices. The interaction terms of the Lagrangian are products of these bilinears. 
In principle, a general effective Lagrangian can be written as a power series 
in the currents ipOrTip and their derivatives. However, it is well known from 
numerous applications of relativistic mean-field models that properties of sym- 
metric and asymmetric nuclear matter, as well as empirical ground state prop- 
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erties of finite nuclei, constrain only the isoscalar-scalar, the isoscalar-vector, 
the isovector-vector, and to a certain extent the isovector-scalar channels. 



The model that we consider in the present work includes the following four- 
fermion interaction vertices: 

isoscalar-scalar: {ipipY 

isoscalar-vector: {tp'~ffj_tlj){tlj'-f^ifj) 

isovector-scalar: {'ipTtp) ■ {ipTtp) 

isovector-vector: {i^if'Juip) ■ {ipf'j'^il') . 

Vectors in isospin space are denoted by arrows. The model is defined by the 
Lagrangian density 

£. = /2free "I" "I" -^der "I" -^erii) (2) 

with the four terms specified as follows: 

Uee = ^{tl^d^-M)i;, (3) 

-\GTvmi'ri,^) ■ {i^f^y) . (4) 

-\ Dvmdu^i.^Wi^r^iJ) + . . . , (5) 

£em = eA^i^T^V^ - -^F.uF^" . (6) 

This Lagrangian (2)-(6) is understood to be formally used in the mean-field 
approximation [19], with fluctuations encoded in density-dependent couplings 
Gi{p) and Di{p), to be specified in detail later. In addition to the free nucleon 
Lagrangian £free and the interaction terms contained in C^f, when applied to 
finite nuclei the model must include the coupling Ccm of the protons to the 
electromagnetic field A^^, and derivative terms contained in £der ^ • One can, of 
course, construct many more derivative terms of higher orders. In this way we 



In the terms duiijjTijj) the derivative is understood to act on both i/j and ip 
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take into account leading effects of finite range interactions that are important 
for a quantitative fine-tuning of nuclear properties. In practice it will turn out 
that a single term Ds{d^ipip){d'^ijj'ip) is already sufficient to account for the 
detailed surface structure of the nuclei considered. 

Our model Lagrangian is formally equivalent to the ones used in the stan- 
dard relativistic mean-field point-coupling models of Refs. [20,21,8,9,10]. The 
underlying dynamics is, however, quite different. The parameters of the in- 
teraction terms in the standard point-coupling models have constant values 
adjusted to reproduce the nuclear matter equation of state and a set of global 
properties of spherical closed-shell nuclei. In order to describe properties of 
finite nuclei on a quantitative level, these models include also some higher 
order interaction terms, such as six-fermion vertices (ipip)^, and eight-fermion 
vertices {ipipY [{il)'^^ip){'4}'^^il))Y . The model that we develop in this work 
includes only second order interaction terms. On the other hand, the coupling 
parameters are not constants, but rather functions of the nucleon density op- 
erator p. Their functional dependence will be determined from finite-density 
QCD sum rules and in-medium chiral perturbation theory. 

Medium dependent vertex functions have also been considered in the frame- 
work of relativistic mean-field meson-exchange models. In the density depen- 
dent relativistic hadron field (DDRH) models the medium dependence of the 
meson-nucleon vertices is expressed as a functional of the baryon field opera- 
tors. The meson-nucleon vertex functions are determined either by mapping 
the nuclear matter Dirac-Brueckner nucleon self-energies in the local density 
approximation [19,22,23], or the parameters of an assumed phenomenological 
density dependence of the meson-nucleon couplings are adjusted to repro- 
duce properties of symmetric and asymmetric nuclear matter and finite nu- 
clei [24,25]. In practical applications of the DDRH models the meson-nucleon 
couplings are assumed to be functions of the baryon density ifj'^'il). In a rela- 
tivistic framework the couplings can also depend on the scalar density 
Nevertheless, expanding in ip'^ip is the natural choice, for several reasons. The 
baryon density is connected to the conserved baryon number, unlike the scalar 
density for which no conservation law exists. The scalar density is a dynamical 
quantity, to be determined self-consistently by the equations of motion, and 
expandable in powers of the Fermi momentum. For the meson-exchange mod- 
els it has been shown that the dependence on baryon density alone provides a 
more direct relation between the self-energies of the density-dependent hadron 
field theory and the Dirac-Brueckner microscopic self-energies [23]. Moreover, 
the pion-exchange contributions to the nucleon self-energy, as calculated us- 
ing in-medium chiral perturbation theory, are directly given as expansions in 
powers of the Fermi momentum. Following these considerations, we express 
the coupling strengths of the interaction terms as functions of the baryon den- 
sity, represented by the operator p = i/j'^iIj in the rest frame of the many-body 
system. 
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2.2 Equation of motion and nucleon self-energies 



The single-nucleon Dirac equation is derived by variation of the Lagrangian 
(2) with respect to ip: 

[^^(,aM_yM)_(M + 5)]V; = 0, (7) 



where 

^/^ = S'^ + r-S^ + E^f (8) 

and 

5 = S5 + r-ST5 + S,s, (9) 
with the nucleon self-energies defined by the following relations 

j:^^ = Gy{^Y^)-eA^^^ , (10) 

^^T = GMi^rri') , (11) 

J:s = Gs{^^J^P)-DsD{^|J^|J) , (12) 

i:TS = GTs{^f^P) , (13) 



-'rS '■ 



dp 



id.f)u,{d^{^ij)) , (14) 



S^f = y ('^(^^)(V^) + ^(^^V^) ■ i^rij) 



dG - - dG 



+ ^{d-^ii^mdu{^i^))^ ■ (15) 

We write pu^ = ^'j'^^ and the four-velocity is defined as (1 — f^)~^/^(l, v) 
where v is the three- velocity vector (f = in the rest-frame of the nuclear 
system). In addition to the isoscalar- vector E'^, isoscalar-scalar E5, iso vector- 
vector and isovector-scalar T,ts self-energies, the density dependence of the 
vertex functions produces the rearrangement contributions E^s and E^. The 
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rearrangement terms result from the variation of the vertex functionals with 
respect to the baryon fields in the density operator p (which coincides with the 
baryon density in the nuclear matter rest-frame). For a model with density de- 
pendent couplings, the inclusion of the rearrangement self-energies is essential 
for energy-momentum conservation and thermodynamical consistency (i.e. for 
the pressure equation derived from the thermodynamic definition and from 
the energy-momentum tensor) [19,24]. 

The density dependence of the coupling strengths which determine the self- 
energies (10-15), will be constrained by in-medium QCD sum rules and chiral 
pion dynamics, to be specified in detail in the following subsections 3.2 and 3.3. 

When applied to nuclear matter or ground-state properties of finite nuclei, the 
point-coupling model is understood to be used in the mean-field approxima- 
tion. The ground state of a nucleus with A nucleons is the antisymmetrized 
product of the lowest occupied single-nucleon self-consistent stationary solu- 
tions of the Dirac equation Eq. (7). The ground state energy is the sum of 
the single-nucleon energies plus a functional of proton and neutron scalar and 
particle densities. 



3 Nuclear matter equation of state 

3. 1 Framework 

In the translationally invariant infinite nuclear matter all terms involving the 
derivative couplings (5) drop out. The spatial components of the four-currents 
vanish, and the densities are calculated by taking expectation values 



where |$) is the nuclear matter ground state. The energy density e and the 
pressure P are derived from the energy-momentum tensor as 



p, = ($|^^|<i>) = p^ + p:, 

p=(<|.|^7Vl$)=/ + p\ 
p,3 = ($|^r3V^|$)=p?-p:, 



(16) 
(17) 
(18) 
(19) 




\ ^ \ 2I2I 2 

P = + - - el,^ + \Gyp' + \GtvpI + \GspI + \GtspI, 
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1 dGs 2 , 1 3 , 1 OGtv 2 , 1 <9Gt5 2 /m \ 

^2^P^P + 2 + 2 ^^^^ + 2 ^^-^ • ^2') 
The particle density is related to the Fermi momentum kj in the usual way, 

- (4 / - ^ ' '''' 



where the index i = p,n refers to protons and neutrons, respectively. The 
corresponding scalar densities are determined by the self-consistency relation 



(27r) 



M* 



Ml 
27r2 



k}E^ 



(M* 



2, k) + E' 



M* 



(23) 



with 



(24) 



and the effective nucleon masses 



= M + GsPs ± Gt5Ps3 • (25) 

The kinetic contributions to the energies of the protons and neutrons in nuclear 
matter are calculated from 

4n=(^ / d'k^WV(m' = \[^E^P' + M:p\], z = p,n. (26) 



Note that, in contrast to the energy density, rearrangement contributions ap- 
pear explicitly in the expression for the pressure. 



3.2 QCD constraints 



We proceed now with a central theme of our work: establishing connections 
between the density- dependent point couplings in (10-15) and constraints from 
QCD. Two key features of low-energy, non perturbative QCD are at the ori- 
gin of this discussion: the presence of a non-trivial vacuum characterized by 
strong condensates and the important role of pionic fluctuations governing 
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the low-energy, long wavelength dynamics according to the rules imposed by 
spontaneously broken chiral symmetry. 



The basic conjecture underlying our model is, consequently, that the nucleon 
isoscalar self-energies arise primarily through changes in the quark condensate 
and in the quark density at finite baryon density, together with chiral (pionic) 
fiuctuations induced by one- and two-pion exchange interactions. In-medium 
QCD sum rules relate the changes of the scalar quark condensate and the 
quark density at finite baryon density, with the scalar and vector self-energies 
of a nucleon in the nuclear medium. In leading order which should be valid 
at densities below and around saturated nuclear matter, the condensate part 
of the scalar self-energy is expressed in terms of the density dependent chiral 
condensate as follows [11,12,13]: 

4°' = I(f.>, - (.-.)ol = ^^P. ■ (27) 

rriu + nid 



The chiral vacuum condensate {qq)o is a measure of spontaneous chiral sym- 
metry breaking in QCD. At a renormalization scale of about 1 GeV (with 
quark masses m„ + ^ 12 MeV) its value [26] is {qq)o — —(240 MeV)^ ~ 
—1.8 fm~^. In order to appreciate the strength of this scalar condensate, note 
that its magnitude is more than a factor of ten larger than the baryon den- 
sity in the bulk of a heavy nucleus. The difference between the vacuum con- 
densate {qq)o and the one at finite density involves the nucleon sigma term, 
(Jat = {N\mqqq\N), to this order. Furthermore, ^ IGeV is a characteristic 
scale (the Borel mass) roughly separating perturbative and non-perturbative 
domains in the QCD sum rule analysis. To the same order in the condensates 
with lowest dimension, the time component of the vector self-energy is 



where the quark baryon density is simply related to that of the nucleons by 

When taken at the same leading order, the QCD sum rule analysis also 
identifies the free nucleon mass at p = according to loffe's formula [27], 
M = —^(qq)o- Together with the Gell-Mann, Oakes, Renner relation 

{niu + rnd){qq)o = —^"if^ (where m^r is the pion mass and = 92.4 MeV is 
the pion decay constant), one finds 

4°^(p)=M*(p)-M = -^p. (29) 

TtJ TV 
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and 



4(m„ + md )M 



-P 



(30) 



Given these self-energies arising from the condensate background, the corre- 
sponding equivalent point-coupling strengths G^y are simply determined by 

= 4%. (31) 



and 



.(0) 
-'V 



P 



(32) 



The second important ingredient is the contribution to the nucleon self-energies 
from chiral fluctuations related to pion-exchange processes (primarily iterated 
one-pion exchange with small corrections from one-pion exchange Fock terms). 
While pion-exchange processes are not explicitly covered by standard Hartree- 
type mean field descriptions, their effects can nevertheless be re-expressed as 
density dependent corrections to the mean fields. The corresponding point- 
coupling strengths, G^gY{p) etc., and their explicit density dependencies are 
calculated using in-medium chiral perturbation theory to a given loop order. 

Combining effects from leading QCD condensates and pionic fluctuations, the 
strength parameters of the isoscalar four-fermion interaction terms in the La- 
grangian (2) are: 

G5y(p) = GgU4:kp)- (33) 



For the isovector channels we assume in this work that only pionic (chiral) 
fluctuations contribute. 

In Ref. [15] the nuclear matter equation of state (EOS) has been calculated us- 
ing in-medium chiral perturbation theory up to three loop order in the energy 
density, expanded in powers of the Fermi momentum kf (modulo functions 
of kf/m.,^). The empirical saturation point, the nuclear matter incompress- 
ibility, and the asymmetry energy at saturation can be well reproduced at 
order 0{k^f) in the chiral expansion with just one single momentum cut-off 
scale A ~ 0.65 GeV which parameterizes short-distance physics. In this ap- 
proach the nuclear matter saturation mechanism is entirely determined by 
in-medium pion-nucleon dynamics: the interplay between attraction from two- 
pion exchange processes and the stabilizing effects of the Pauli principle. The 
implicit assumption made in the analysis of Ref. [15] is the following: there is 
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no net contribution to the energy density around saturation from the nucleon 
isoscalar self-energies, S^'* and T,y\ that arise through changes in the quark 
condensate and the quark density at finite baryon density. At first sight, this 
assumption appears to be unreahstic, given that in-medium QCD sum rules 
suggest individually large scalar and vector self-energies of about 300 - 400 
MeV in magnitude from the condensate background (see Eqs. (29)- (30)). On 
the other hand, the ratio 

= ^ Pi (34) 

4(m„ + md) p 



is indeed consistent with —1, for typical values of the nucleon sigma term 
(Tn and the current quark masses rriu^d, and around nuclear matter saturation 
density where ps — p (as an example, take aN — 50 MeV [28] and mu+rrid ~ 12 
MeV [26] at a scale of IGeV). One should note of course that the QCD sum rule 
constraints implied by (29)- (30) and by the ratio (34) are not very accurate 
at a quantitative level. The leading-order loffe formula on which Eq. (29) 
relies has corrections from condensates of higher dimension which are not well 
under control. The estimated error in the ratio T,^ /^P ~ -1 is about 20%, 
given the uncertainties in the values of and m„ + m^. Nevertheless, the 
constraints implied by Eq. (34) give important hints for further orientation. 



3.3 Equation of state based on in-medium Chiral Perturbation Theory 



In this work the nuclear matter equation of state will be calculated in three 
steps. In the first approximation we follow the implicit assumption made in 
Ref. [15] that E^^ = —Tly^ in nuclear matter, and neglect the contribution 
of the condensate background self-energies to the nuclear matter equation of 
state (EOS). The density dependence of the strength parameters is determined 
by equating the isoscalar-scalar, the isoscalar-vector, the isovect or- scalar, and 
the isovector-vector self-energies Eqs. (10) - (15) in the single-nucleon Dirac 
equation (7) with those calculated using in-medium chiral perturbation theory: 



G^s^Ps = ^T\k,^p)^ (35) 

GWp + SW = S^«PT(A;„p), (36) 

41p.3 = Sgr^(A;/,p), (37) 

GWp3 = S?«PT(fc;,p), (38) 



with 
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where the superscripts vr indicates that the density dependence of the cou- 
phngs originates from pionic (chiral) fluctuations. The EOS of symmetric and 
asymmetric nuclear matter calculated in CHPT give, via the Hugenholtz-van 
Hove theorem, the sums U(T){kf,kf) = S|^|'^(/cj, p) + H'fj^^^kf, p) of the 
scalar and vector nucleon self-energies in the isoscalar and isovector channels 
at the Fermi surface p = /c/ up to two- loop order, generated by chiral one- 
and two-pion exchange [16]. The difference T,'^r^^^{kf, p) — E^y'^(/c/, p) is cal- 
culated from the same pion-exchange diagrams via charge conjugation of the 
single particle potential in nuclear matter, as explained in Ref. [16]. Following 
a procedure similar to the determination of the nucleon-meson vertices of rel- 
ativistic mean- field models from Dirac-Brueckner calculations [23], we neglect 
the momentum dependence of SpEjgy (p, p) and take their values at the Fermi 
surface p = kf. A polynomial fit up to order k^ is performed, and all four 
CHPT self-energies have the same functional form. 



E(/c;,A) 



C30 + C31A C32A^ + C3L In 
C50 + C5L In 



+^40 ^ + 



M2 



47r/^A 



1,5 



(40) 



where the dimensionless parameter A is related to the momentum cut-off scale 
by A = 27r/7rA, and f.^ = 92.4 MeV denotes the pion decay constant. In Table 1 
we display the coefficients of the expansion (40) of the CHPT self-energies for 
all four channels. The cut-off parameter A = 646 MeV that has been adjusted 
to the nuclear matter saturation in Ref. [15], corresponds to A ~ 1.113. In 
order to determine the density dependence of the couplings from Eqs. (35) - 
(38), the CHPT self-energies are re-expressed in terms of the baryon density 
p = 2H/37r2: 



Eg"PT(A:^,p) = (c,i + c,2P^ + c,3p3)p , (41) 

E^«PT(/cj, p) = (c„i + c^2P^ + c„3P^)p , (42) 

^rT'^ikf, P) = (Ctsl + Cts2p^ + Q,3P^)P3 , (43) 

^Tv'^ikf, P) = {Ctvl + Ctv2p^ + Cto3pi)p3 . (44) 

The resulting expressions for the density dependent couplings of the pionic 
fluctuation terms are ^ 



G^S^ = Csl + C,2P^ + C,3p^ , (45) 



2 small differences between ps and p at nuclear matter densities are neglected here 
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j{QCy2 - Cs2 - S\cts2 + Ctos)) . (46) 
i(3c„3 - Cs3 - 5^(Qs3 + Ctos)) 

(47) 
(48) 

where 6 = {pn — Pp)/{pn + Pp)- For A ~ 1.113, the coefficients of the expan- 
sion of the CHPT self-energies in powers of the baryon density Eqs. (41) - 
(44) are given in Table 2. In Fig. 1 the resulting equation of state of isospin 
symmetric nuclear matter, calculated from Eqs. (20) and (21), is compared 
with the CHPT nuclear matter EOS of Ref. [15]. In Table 3 we compare the 
corresponding nuclear matter properties at the saturation point: the binding 
energy per particle, the saturation density, the compression modulus, and the 
asymmetry energy. The point-coupling model, with the density dependence of 
the strength parameters of the second order interaction terms determined by 
pionic (chiral) fluctuations, nicely reproduces the nuclear matter EOS calcu- 
lated using in-medium CHPT. Small differences arise mainly because in the 
mapping of the CHPT nucleon self-energies on the self-energies of the point 
coupling model in Eqs. (35) - (38), the momentum dependence of the former 
has been frozen to the values at the Fermi momentum. This is a well known 
problem that arises also in the determination of the meson-nucleon in-medium 
vertices from Dirac-Brueckner calculations of nucleon self-energies in nuclear 
matter [23]. For comparison in Fig. 1 we also include the microscopic many- 
body nuclear matter equation of state of Friedman and Pandharipande [29]. 

The importance of the rearrangement self-energies implemented through Eqs. 
(36) and (46) for a model with density dependent couplings is illustrated in 
Fig. 2, where we compare the EOS of isospin symmetric nuclear matter cal- 
culated with and without the rearrangement contributions in the Eq. (36) 
that determines the density dependence of the vertex of the isoscalar-vector 
interaction term. The difference is indeed very large. Without the rearrange- 
ment terms in Eq. (36) a saturation density of only 0.1 fm~^ is obtained, 
and the binding energy per particle is —6.85 MeV. Also the nuclear matter 
incompressibility is too low: 124 MeV. Obviously, only with inclusion of the 
rearrangement contributions in the isoscalar-vector self-energy it is possible 
to reproduce the nuclear matter EOS calculated in CHPT in an equivalent 
approach with density dependent point couplings. 

At very low-densities, the CHPT equation of state (or any mean-field ap- 
proach) cannot reproduce realistic calculations of E/A, for obvious reasons. 
The low density limit of the energy per particle is not given by a free Fermi gas 
of nucleons, but by a gas of clusters (deuterons, etc. . .). As a consequence, the 
low-density limit of E/A approaches a constant as kf 0. The effect of this 



G 



(tt) 
V 



C„i + Cy2p-i + C^sP'^ 



Cv3 



Gj^Y = Ctvl + Ctv2P'* + CtvSp^ 
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shift is marginal, however, for the bulk of nuclear matter around saturation 
density. 



3.4 Constraints from QCD Sum Rules: leading orders 



The next step in our calculation of the nuclear matter EOS includes the con- 
tributions of the condensate background self-energies in the isoscalar-scalar 
and isoscalar- vector channels. In Ref. [17] we have shown that, even though in 
a first approximation the condensate potentials do not play a role in the satu- 
ration mechanism, they are indeed essential for the description of ground-state 
properties of finite nuclei. 

The nuclear matter EOS is calculated from Eqs. (20) and (21). In this case, 
the strength parameters of the isoscalar-scalar and isoscalar-vector interaction 
terms have the form given in Eq. (33), i.e. they include contributions from 
condensate background fields. For the isovector channel, on the other hand, 
the assumption of the model is that only pionic (chiral) fluctuations contribute 
to the nucleon self-energies. In leading order of the in-medium QCD sum rules 
the nucleon isoscalar self-energies that arise through changes in the quark 
condensate and the quark density are linear functions of the corresponding 
scalar and baryon densities (see Eqs. (31) and (32)). To this order, therefore, 
the corresponding couplings should be density independent. Eq. (29) 

implies (identifying M with the free nucleon mass): 

T.f ^ -350 MeV ^ . (49) 

^ 50 MeV po 



Evidently, the leading-order in-medium change of the chiral condensate is a 
source of a strong, attractive scalar field which acts on the nucleon in such a 
way as to reduce its mass in nuclear matter by more than 1/3 of its vacuum 
value. Using (31) one estimates: 

4°) ^ -11 fm^ at p,~po = 0.16fm-^ (50) 



A correspondingly smaller value of Gg results if only a fraction of the nucleon 
mass M is associated with the chiral condensate, leaving room for non-leading 
contributions from higher dimensional condensates. Since the QCD sum rules 
suggest that the ratio of the condensate scalar and vector self-energies is close 
to —1, one expects roughly Gy^ ~ ~G^P ■ 

In our model both and Gy^ are parameters to be adjusted to "empirical" 
saturation properties of nuclear matter, as well as to ground-state properties 



16 



of finite nuclei. Moreover, since the nuclear matter saturation mechanism is 
not any more exclusively determined by pionic (chiral) fluctuations, the mo- 
mentum cut-off scale A (or the dimensionless A) is an additional parameter in 
our calculation of the nuclear matter EOS. At this preliminary stage, the three 
parameters: G^p , Gy^ and A are adjusted to reproduce "empirical" nuclear 
matter properties: E/A = -16 MeV (5%), po = 0.153 fm'^ (10%), Kq = 250 
MeV (10%), and = 33 MeV (10%). The values in parentheses correspond 
to the error bars used in the fitting procedure. 

The result of the corresponding least-squares fit is shown in Fig. 3 (EOS of 
isospin symmetric nuclear matter in comparison with the one calculated exclu- 
sively from chiral pion-exchange) and in the first row of Table 4 for G^ = -7 
fm^, Gy^ = 7 fm^, and A = 685 MeV. The EOS so obtained is not yet satisfac- 
tory. The asymmetry energy at saturation is too high, and the values of G^P 
and are smaller than the leading order QCD sum rule estimate Eq. (50). 
This results in the relatively large Dirac effective mass M* /M = 0.75. A high 
Dirac mass in turn indicates that the effective single-nucleon spin-orbit poten- 
tial is too weak to reproduce the empirical energy spacings between spin-orbit 
partner states in finite nuclei. Stronger background scalar and vector fields are 
required in order to drive the large spin-orbit splittings in finite nuclei. 

Notice however that, even though G^g^ and have been varied indepen- 
dently, the minimization procedure prefers to balance the contributions from 
the corresponding large scalar and vector self-energies so that their sum tends 
to vanish. There is already enough binding from pionic fluctuations alone, 
and therefore E^'* = — Sy"* represents a very good approximation for the con- 
densate potentials. This remarkable feature prevails as we now proceed with 
fine-tuning improvements. 



3. 5 Corrections of higher order 



Up to this point the nucleon self-energies are evaluated at two-loop level in- 
cluding terms of order /cj. At order /cj (oc p^), several additional effects appear. 
First, condensates of higher dimension enter, such as (^EggEg), {q^q)"^ and the 
gluon condensate, the detailed in-medium dependence of which is difficult 
to estimate. Secondly, four-loop CHPT contributions to the energy density 
(three-loop in the self energies) introduce genuine 3-body interactions. All 
those effects combine to produce the 0{k^) correction in the energy per parti- 
cle. Our chiral counting conjecture is that these terms should be subleading in 
the /c/-expansion. In estimating those contributions, we will take a pragmatic 
point of view. Given the density-independent leading condensate terms G^'^^ 
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and the CHPT couplings G^'^\p) evaluated to O(fcj) oc pa, we generalize 
G{p) = + (p) + 5G^^^ (p) , (51) 

adding terms ^G*-^^ = (^'-^V ^^t the constant g^^^ be determined by a least 
squares fit to properties of nuclear matter and finite nuclei. At first sight, the 
corrections 5G'g\ 5Gy \ etc .. . would then introduce several additional pa- 
rameters, undermining our startegy to keep the number of freely adjustable 
fine-tuning constants as small as possible. For the scalar self-energy a primary 
uncertainty arises from the four-quark condensate and its assumed factoriza- 
tion into the form (qq)"^, as discussed in detail in Ref. [13]. While the in-medium 
values of condensates such as {q^qq^q) are not well determined, it neverthe- 
less appears that only a weak or moderate density dependence is consistent 
with known nuclear phenomenology. We do not have to be much concerned 
about this issue since our approach explicitly includes scalar vr — tt fluctuations 
(up to order fcj in the energy density), which should account for a large part 
of the four quark condensate effects. Remaining short-distance effects can be 
absorbed in 6G'^^\ 

As a first option, assume that the only non-vanishing ^G*-^-* is 5Gy^ = gy^ p, 
and determine the constant (y'y'*. For G^g^ = —12 fm^, Gy^ = 11 fm^, gy^ = 
—3.9 fm^ and A = 600 MeV, determined in a least- squares fit to the "empirical" 
nuclear matter input, the EOS resulting is displayed in Fig. 3 and the nuclear 
matter properties at saturation are listed in the second row of Table 4. This 
nuclear matter EOS is actually quite satisfactory, with realistic values of the 
binding energy, saturation density and a low effective Dirac mass. The only 
exception is a relatively large nuclear matter incompressibility. Calculations 
of the excitation energies of isoscalar giant monopole resonances in spherical 
nuclei in the time-dependent relativistic mean-field framework [30] , and in the 
relativistic random-phase approximation [31], suggest that the nuclear matter 
incompressibility modulus should be in the range Kq ^ 250—270 MeV. A value 
of Kq in this range, however, can only be obtained by including the next higher 
order non-linear self-interaction term in the expansion of Sy or Eg, i.e. a term 
proportional to p^. This feature is also well known in standard meson-exchange 
and point-coupling relativistic mean-field models. In the present analysis we do 
not aim for a detailed description of nuclear properties that depend crucially 
on the incompressibility, and we also wish to keep the number of adjustable 
parameters at minimum, so the density dependence of Gy is parameterized as 

Gy{p) = GP+G^^\p)+g'^^p. (52) 

The resulting couplings of the condensate background fields {G^ = —12 fm^ 
and Gy^ = 11 fm^) are remarkably close to the prediction (50) of the leading 
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order in-medium QCD sum rules. The scalar and the vector "condensate" self- 
energies S^"* = Gg Ps and Hy^ = Gy^p follow the expectation t!^ ^ — Sy"* 
to within 5% at saturation density, even without this condition being pre- 
imposed. The large isoscalar condensate background self-energies in turn lead 
to a relatively low effective Dirac mass, crucial for the empirical spin-orbit 
splittings in finite nuclei. Finally, the cut-off A = 600 MeV differs by less than 
10 % from the value (646 MeV) obtained when the nuclear matter EOS results 
exclusively from one- and two-pion exchange between nucleons [15]. 

The SGy^ = gy^p term acts like a three-body force in the energy density. Its 
effect is relatively small: at saturation density, 5Gy^ /Gy^ ~ 0.05. Splitting 
5G^^^ between scalar and vector parts (e.g. by choosing 5Gy^ ~ 5G^P ~ 
—2 fm~^) would induce only marginal changes. 



3.6 Asymmetric nuclear matter 



The energy per particle of asymmetric nuclear matter can be expanded about 
the equilibrium density psat in a Taylor series in p and 5 [32] 

E(p, 5) = E{p, 0) + S2{p)S^ + S,{p)6' + ■■■ (53) 
where 5 = (p„ - Pp)/ (p„ + Pp), and 



E{p, 0) = -a, + (p - psat)' + ... , (54) 



S2{p) = a4 + ^ (p - Psat) + (P - Psat)^ + " " " . (55) 

Psat J-OPsat 



The empirical value at saturation density is S'2(psat) = 04 = 30 ±4 MeV. The 
parameter Pq defines the linear density dependence of the symmetry energy, 
and AKq is the correction to the incompressibility. The contribution of the 
term S'4(p)a^ in (53) is very small in ordinary nuclei and it will be neglected in 
the present work. In a recent analysis of neutron radii in non-relativistic and 
covariant mean-field models [33], Furnstahl has studied the linear correlation 
between the neutron skin and the asymmetry energy. In particular, he has 
shown that there is a very strong linear correlation between the neutron skin 
thickness in 208p]^ ^^^^ ^^le individual parameters that determine the asjmime- 
try energy 5*2 (p): 04, po and AKq. The empirical value of the difference between 
neutron and proton rms radii r„ — Vp in ^osp]^ (0.20 ± 0.04 fm from proton 
scattering data [34], and 0.19 ± 0.09 fm from the alpha scattering excitation 
of the isovector giant dipole resonance [35]) places the following constraints 
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on the values of the parameters of the symmetry energy: 04 ~ 30 — 34 MeV, 
2 MeV/fm^ < po < A MeV/fm^, and -200 MeV < AKq < -50 MeV. 

In Fig. 4 we plot the asymmetry energy 5*2 (p) as function of the baryon density 
for the three steps of our calculation of the nuclear matter EOS: first, using 
only chiral one- and two-pion exchange calculated in CHPT, next with the in- 
clusion of leading order isoscalar condensate background nucleon self-energies 
that arise through changes in the quark condensate and the quark density 
at finite baryon density, and finally with the non-linear contribution 5Gy^ 
to the vector condensate nucleon self-energy. The three curves can be com- 
pared, for instance, with the asymmetry energy obtained in the non-relativistic 
Brueckner-Hartree-Fock calculation of Ref. [36], or with the asymmetry en- 
ergy of the relativistic meson-exchange model with a phenomenological density 
dependence of the meson-nucleon couplings in Ref. [25] . 

In Fig. 5 we display the energy per particle of neutron matter as a function of 
the neutron density. The results obtained from chiral one- and two-pion ex- 
change between nucleons, by adding the isoscalar condensate background nu- 
cleon self-energies linear in the corresponding densities, and finally by includ- 
ing the non-linear "3-body" contribution to the isoscalar vector self-energy, 
are shown in comparison with the microscopic many-body neutron matter 
equation of state of Friedman and Pandharipande [29] . 

These comparisons show a qualitatively correct trend in both the asymmetry 
energy and the neutron matter equation of state, but refinements of the density 
dependent isovector couphngs are obviously still necessary. 

Of course, one cannot expect to reproduce the behaviour at very low neutron 
density in a perturbative approach. The very large neutron-neutron scattering 
length requires a partial resummation which effectively renormalizes the slope 
of the kj-teim in E^/N by roughly a factor of one-half [33,37]. 



3.7 Intermediate summary and discussions 



While the quantitative details for the terms contributing to the scalar and 
vector self-energies can be reconstructed from the previous sections using Ta- 
bles 1 and 2, it is instructive to observe systematic trends in these self-energies 
when expanded in powers of kf (or equivalently, in powers of ps). We focus 
on symmetric nuclear matter and write 

Es,v = + ^s,v + S^s,v , (56) 
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combining the leading condensate terms S^y, the chiral (pionic) terms 
calculated to order kj, and possible corrections of higher order summarized 
in STisy- We introduce the saturation density po = 0.16 fm~'^ as a convenient 
scale and write the condensate terms as 



eS^^(p)~+0.35 GeV 



(P) 



-0.35 GeV 







11 fm^ 


\Po 


1-^(0)1 
I'^V \ 


i- 


11 w 





(57) 
(58) 



We note that with \G 



chosen less than 5% larger than Gy\ this just com- 



pensates for the difference between the baryon density p and the scalar den- 
sity ps < p such that S^^ ~ — S^'* results at p = Pq. The "best-fit" values 

+ 11 fm^ are thus perfectly consistent with QCD 



G 



(0) 



-12 fm^ and gJ?^ 



sum rule expectations: a non-trivial result. The pionic fluctuations from one- 
and two-pion exchange processes have the following approximate pattern: 



'S,V 



(p) ~ -75 MeV ( — 




(59) 



where di varies between —0.61 (scalar) and —0.65 (vector), and d2 — —0.17 
for both cases. Finally, the higher order corrections are summarized as 




SEs + 5Ev - -20 MeV — , (60) 



where the detailed decomposition into scalar and vector part turns out not to 
be relevant. 

The leading attraction at 0{kj) in Ti^gy depends on the cut-off scale A = 
27r/7rA which distinguishes "active" two-pion exchange dynamics at long and 
intermediate distances from short- distance (high momentum) dynamics. De- 
tails at short distance scales (related to the intrinsic structure of the nucle- 
ons and their short-range interactions) are not resolved as long as the Fermi 
momentum kf is much smaller than ATifj^. The unresolved short-distance in- 
formation can then be translated into an equivalent (density-independent) 
four-point vertex which generates self-energy terms linear in the density at 
the mean-field level. 
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3.8 Comparison with the Dirac-Brueckner G matrix 



At this point it is instructive to compare the density- dependence of our self- 
energies T,sy{p) with that of the scalar and the vector self-energies resulting 
from full Dirac-Brueckner calculations based on realistic NN-interactions. We 
refer here explicitly to a calculation, reported in Ref. [38], which uses the Bonn 
A potential. With this potential, the nuclear matter saturation density comes 
out somewhat higher {psat — 0.185 fm~^) than the one of our best fit. Simu- 
lating this higher saturation density in our approach requires weakening the 
condensate mean fields S^y by about 25 %, but with no changes in the pionic 

terms Sg^y. After this readjustment, the difference in the /c/-dependences be- 
tween our T,s and and those resulting from the Dirac-Bruckner calculation 
is less than 10 % over the entire range of densities from 0.5po to 2.5po (see 
Fig. 6). 



This is a remarkable observation: it appears that in-medium chiral perturba- 
tion theory at two-loop order, with a cut-off scale A 2^ 0.6 GeV that converts 
unresolved short-distance dynamics at momenta beyond A effectively into con- 
tact terms, generates quantitatively similar in-medium nucleon self energies 
as a full Dirac-Brueckner calculation, when requiring that both approaches 
reproduce the same nuclear matter saturation point. The reasoning behind 
this observation can presumably be traced back to the separation of scales at 
work, a key element of chiral effective field theory on which our approach is 
based. One- and two-pion exchange in-medium dynamics at momentum scales 
comparable to kf are treated explicitly. Our in- medium CHPT approach in- 
cludes all terms (ladders and others) to three-loop order in the energy density, 
a prominent one being iterated one-pion exchange. The non-trivial kf depen- 
dence in the self-energies (beyond "trivial" order fcj) reflects the action of 
the Pauli principle in these processes. Such features are also present in lead- 
ing orders of the Brueckner ladder summation which includes, for example, 
Pauli blocking effects on iterated one-pion exchange. These terms produce the 
characteristic /c|-dependence of the energy per particle. On the other hand, 
the iteration of short- distance interactions to all orders in the Brueckner lad- 
der, involves intermediate momenta much larger than kj and generates similar 
self-energy pieces (proportional to density p) as a contact interaction reflecting 
the high-momentum scale A. Note that this contact interaction must not be 
iterated as it already represents the full short-distance T-matrix information. 



The reasoning here has a close correspondence to related studies in Ref. [39]. 
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4 Finite nuclei 



In this section the relativistic nuclear point-coupling model, constrained by 
in-medium QCD sum rules and chiral perturbation theory, is applied in cal- 
culations of ground state properties of finite nuclei. In the mean-field ap- 
proximation the ground state of a nucleus with A nucleons is represented by 
the antisymmetrized product of the lowest occupied single-nucleon stationary 
solutions of the Dirac equation (7). The calculation is self-consistent in the 
sense that the nucleon self-energies are functions of nucleon densities and cur- 
rents, calculated from the solutions of Eq. (7). The mean-field approach to 
nuclear structure represents an approximate implementation of Kohn-Sham 
density functional theory (DFT) [40,41,42], which is successfully employed 
in the treatment of the quantum many-body problem in atomic, molecular 
and condensed matter physics. At the basis of the DFT approach are en- 
ergy functionals of the ground-state density. In relativistic mean-field models, 
these become functionals of the ground-state scalar density and of the baryon 
current. The scalar and vector self-energies play the role of local relativistic 
Kohn-Sham potentials [1,3]. The mean- field models approximate the exact en- 
ergy functional, which includes all higher-order correlations, with powers and 
gradients of densities, with the truncation determined by power counting [3]. 

In this work we only consider spherical even-even nuclei. Because of time- 
reversal invariance the spatial components of the four-currents vanish, and 
the nucleon self-energies (reduced to the time component T,^ of the vector 
terms E'^) read: 



^^ = Gvp-Vc, (61) 

S0=GwP3, (62) 

J:s = GsPs + DsAps, (63) 

St5 = GtsPss , (64) 

S.5 = ^(Vp.)-(Vp), (65) 



(66) 



where Vc is the Coulomb potential. In the mean-field approximation the 
isoscalar-scalar, isoscalar- vector, isovector-scalar and isovector-vector ground 
state densities are calculated from the Dirac wave functions ipk of the occupied 
positive-energy orbits as: 
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A 

Ps = Yl ^k^k , (67) 



fc=l 
A 



p = J2mk, (68) 

k=l 

A 

Ps3 = '^k'T^lpk , (69) 



fc=l 
A 



P3 = T.4r3^k, (70) 



k=l 



respectively. The expression for the ground state energy of a nucleus with A 
nucleons reads: 



E = j2ek-^Jd^ GspI + GtspI, + Gyp' + +GtvpI + p'Vc 

k=l 

dGs 2 , 9Gv 3 , OGtv 2 , ^^G'ts 2 , n a 
^^P^P^^P +^P^P+^PssP + DsPsAp, 

^-(Vp.) ■ (Vp.)p + 2%^(Vp.) ■ (Vp)pJ , (71) 



ap ap 

where denotes the single-nucleon energies. In addition, for open shell nuclei 
pairing correlations are included in the simple BCS approximation [43]. After 
obtaining the solution of the self-consistent Dirac equation, the microscopic 
estimate for the center-of-mass correction 

Earn = , (72) 

2AM ' ^ ' 



is subtracted from the total binding energy. Here P^^ is the squared total 
momentum of the nucleus with A nucleons. 

The details of the calculated properties of finite nuclei will, of course, depend 
on the fine-tuning of the coupling parameters. The density-dependent cou- 
plings G5, Gy, Gts and Gtv are constrained by in-medium QCD sum rules 
and chiral perturbation calculations of one- and two-pion exchange diagrams, 
as described in the previous section. For finite nuclei we must also determine 
the coupling parameter of the derivative term: Ds- Dimensional considerations 
suggest the following ansatz 

Os(p) = ^ . (73) 

where M. is a. characteristic mass scale for a given spin-isospin channel. There 
is no deeper reason, however, for the parameters of the derivative terms to have 
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the same density dependence as the couphng parameters of the four-fermion 
interaction terms. The simplest option, followed in this work, is to introduce 
only Ds and treat it as a density-independent adjustable parameter. In this 
case the remaining rearrangement contribution to the vector-self energy be- 
comes, of course, much simplified. As it has been emphasized by Serot and 
Furnstahl [6], the empirical data set of bulk and single-particle properties of 
finite nuclei can only constrain six or seven parameters in the general expan- 
sion of the effective Lagrangian in powers of the fields and their derivatives. In 
particular, only one parameter of the derivative terms can be determined by 
the binding energies and radii of spherical nuclei. We adjust the single remain- 
ing surface parameter Dg of the isoscalar-scalar derivative term to properties 
of light and medium-heavy N ^ Z nuclei. This approximation, which was 
first used by Serot and Walecka in Ref. [44], implies that the isoscalar- vector, 
the isovector-scalar and the isovector-vector interactions are considered to be 
purely contact interactions (no gradient terms). 

Before we present the results of full model calculations for finite nuclei, it is 
instructive to consider separately the contributions of chiral pion dynamics and 
condensate background self-energies to properties of finite nuclei. In the first 
step we have calculated the ground states of ^^O and '^^Ca by using the coupling 
parameters determined by the nuclear matter EOS of Ref. [15] (see Fig. 1): 

Gsip) = G^ip), Gy{p) = G^^\p), Grsip) = Gtv{p) = G^tUp). 

A = 646.3 MeV, while the couplings to the condensate background fields are 
set to zero. In this case the nuclear dynamics is completely determined by 
chiral (pionic) fluctuations. The calculated total binding energies are already 
within 5— 8 % of the experimental values, but the radii of the two nuclei are too 
small (by about 0.2 fm). This is because the spin-orbit partners (lp3/2, lpi/2) 
and (1^5/2, 1*^3/2) are practically degenerate. In Fig. 7 we display the calculated 
neutron and proton single-particle levels of ^^O and ^°Ca. The energies of the 
degenerate doublets are close to the empirical positions of the centroids of the 
spin-orbit partner levels, and even the calculated energies of the s-states are 
realistic. This is a very interesting result. The chiral pion dynamics provides 
the saturation mechanism and binding of nuclear matter, but not the strong 
spin-orbit force. The inclusion of the isoscalar-scalar derivative term has some 
effect on the calculated radii but it cannot remove the degeneracy of the spin- 
orbit doublets. 

The spin-orbit potential plays a central role in nuclear structure. It is at the 
basis of the nuclear shell model, and its inclusion is essential in order to repro- 
duce the experimentally established magic numbers. In non-relativistic models 
based on the mean-field approximation, the spin-orbit potential is included in 
a purely phenomenological way, introducing the strength of the spin-orbit 
interaction as an additional parameter. Its value is usually adjusted to the 
experimental spin-orbit splittings in spherical nuclei, for example ^^O. On the 
other hand, in the relativistic description of the nuclear many-body prob- 
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lem, the spin-orbit interaction arises naturally from the scalar-vector Lorentz 
structure of the effective Lagrangian, and relativistic models reproduce the 
empirical spin-orbit splittings. In the first order approximation, and assuming 
spherical symmetry, the spin-orbit term of the effective single-nucleon poten- 
tial can be written as 

where the large spin-orbit potential Vis arises from the difference of the vector 
and scalar potentials, — S* ~ 0.7 GeV [2,45]. Explicitly, 

Vis = ^iV'-S). (75) 



where Mg// is an effective mass specified as [2] 

Meff = M-^iV'-S). (76) 



The isoscalar nucleon self-energies generated by pion exchange are not suffi- 
cently large to produce the empirical effective spin-orbit potential. The degen- 
eracy of spin-orbit doublets is removed by including the isoscalar condensate 
background nucleon self-energies that arise through changes in the quark con- 
densate and the quark density at finite baryon density. 

In Table 5 we display the binding energies per nucleon and charge radii of light 
and medium-heavy nuclei, calculated in the relativistic point-coupling model 
constrained by in-medium QCD sum rules and chiral perturbation theory, in 
comparison with experimental values. In addition to the four parameters deter- 
mined by the nuclear matter equation of state, we have adjusted the isoscalar- 
scalar derivative term in the calculation of ^^O and ^°Ca: Ds = —0.713 fm^. 
This value of Ds is very close to the ones used in the effective interactions of 
the standard relativistic point-coupling model of Ref. [10]. It is also consistent 
with the "natural" order of magnitude expected from Ds ~ Gs/A."^- The re- 
sulting agreement between the calculated and empirical binding energies and 
charge radii is indeed very good. 

The effect of including the isoscalar condensate background nucleon self-energies 
is illustrated in Figs. 8 and 9, where we display the neutron and proton single- 
particle energies in ^^O and ^°Ca, respectively, in comparison with the cor- 
responding experimental levels. For both nuclei the calculated levels are in 
excellent agreement with the empirical single-nucleon levels in the vicinity of 
the Fermi surface. In particular, with the inclusion of the isoscalar condensate 
self-energies, the model reproduces the empirical energy differences between 
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spin-orbit partner states. This is a very important result and it shows that, 
while nuclear binding and saturation are partially generated by chiral (two- 
pion exchange) fluctuations in our approach, strong scalar and vector flelds of 
equal magnitude and opposite sign, induced by changes of the QCD vacuum 
in the presence of baryonic matter, generate the large effective spin-orbit po- 
tential in flnite nuclei. Not surprisingly, the ^°Ca spectrum is reminiscent of 
an underlying pseudo-spin symmetry [46]. 

The charge density distributions and charge form factors of ^^O and ^°Ca are 
displayed in Figs. 10 and 11, respectively. The results obtained in the rela- 
tivistic point-coupling model, constrained by in-medium QCD sum rules and 
chiral perturbation theory, are compared with the experimental charge den- 
sity proflles and form factors [47] . The theoretical charge density distributions 
are obtained by folding the point proton density distribution with the electric 
nucleon form factor of dipole form. The corresponding form factors are plotted 
as functions of the momentum q [4] 



where pchio) is the Fourier transform of the charge density and the second term 
is a center-of-mass correction. Higher-order effects in /Pcm negligible. We 
notice the excellent agreement between the calculated and experimental charge 
form factors for momenta |g| < 2 fm~^. When the calculated charge densities 
are compared to results obtained with standard relativistic mean-fleld effective 
interactions, the present results display less pronounced shell effects, in better 
agreement with empirical data. 

In Fig. 12 we compare the single-nucleon spectra of ^^Ni, calculated in our ap- 
proach, with the results of a relativistic mean-fleld calculation using NL3 [48], 
probably the best non-linear meson-exchange effective interaction. The spectra 
are in remarkable agreement. All these results demonstrate that in the present 
approach, based on QCD sum rules and in-medium chiral perturbation theory, 
and with a small number of model parameters determined directly by these 
constraints, it is possible to describe symmetric and asymmetric nuclear mat- 
ter, as well as properties of flnite nuclei, at a quantitative level comparable 
with phenomenological relativistic mean-field models. 



5 Summary and conclusions 

The goals of nuclear structure theory have evolved, especially in the last 
decade, from the macroscopic and microscopic descriptions of structure phe- 




(77) 
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nomena in stable nuclei, towards more exotic nuclei far from the valley of (3- 
stability, and towards the investigation of the bridge between low-energy QCD 
and the dynamics of nuclei. New experimental information on exotic nuclei 
with extreme isospins present serious challenges for traditional nuclear struc- 
ture models. It is likely that radically improved models, or even completely 
new theoretical strategies and approaches, will be required as one proceeds to 
these unexplored region with radioactive nuclear beams. 

Models based on effective field theories and density functional theory provide 
a systematic theoretical framework for studies of nuclear dynamics. By using 
the important active degrees of freedom at low energies, and by expanding 
in small parameters determined by the relevant scales, effective Lagrangians 
are constructed that provide a microscopically consistent, and yet simple and 
economical approach to the nuclear many-body problem. 

In this work we have derived, following Ref. [17], a relativistic nuclear point- 
coupling model with density-dependent contact interactions, which emphasizes 
the connection between nuclear dynamics and two key features of low-energy, 
non-perturbative QCD: the presence of a non-trivial vacuum characterized by 
strong condensates and the important role of pionic fluctuations governing 
the low-energy, long wavelength dynamics according to the rules imposed by 
spontaneously broken chiral symmetry. 

Investigations based on QCD sum rules have shown how large isoscalar scalar 
and vector nucleon self-energies, that characterize the nuclear ground state, 
arise in finite-density QCD through changes in the quark condensate and the 
quark density at finite baryon density. Chiral (pionic) fluctuations superim- 
posed on the condensate background fields play a prominent role for nuclear 
binding and saturation. These pionic fluctuations are calculated according to 
the rules of in-medium chiral perturbation theory. 

Based on these observations, we have formulated a density-dependent point- 
coupling model for nuclear matter and finite nuclei. The coupling parameters 
are functions of the nucleon density operator p. Their functional dependence 
is determined from finite density QCD sum rules and in-medium chiral per- 
turbation theory. The strength parameters of the second order isoscalar inter- 
action terms in the Lagrangian include contributions from condensate back- 
ground fields and pionic (chiral) fluctuations. For the isovector channel it is 
(so far) assumed that only the pionic fluctuations contribute to the nucleon 
self-energies. 

The resulting nuclear point-coupling model, formally used in the mean-field 
approximation, but encoding effects beyond mean-field through the density 
dependent couplings, has been employed in an analysis of the equations of 
state for symmetric and asymmetric nuclear matter, and of bulk and single- 
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nucleon properties of light and medium N ^ Z nuclei. In comparison with 
purely phenomenological relativistic and non-relativistic mean-field models of 
nuclear structure, the built-in QCD constraints and the explicit treatment 
of one- and two-pion exchange significantly reduce the freedom in adjusting 
parameters and functional forms of density-dependent vertex functions. 

A very good equation of state of symmetric nuclear matter is obtained already 
by equating the isoscalar-scalar, the isoscalar-vector, the isovect or- scalar, and 
the isovector-vector nucleon self-energies in the single-nucleon Dirac equation 
with those calculated using in-medium chiral perturbation theory up to three 
loop order in the energy density. The contributions to the energy density 
originate exclusively from one- and two-pion exchange between nucleons, and 
they are ordered in powers of the Fermi momentum kf (modulo functions of 
kf/rriTj). The empirical saturation point, the nuclear matter incompressibility, 
and the asymmetry energy at saturation can be well reproduced at order 0{k^) 
in the chiral expansion with just one single momentum cut-off scale A ~ 0.65 
GeV which parameterizes unresolved short- distance dynamics. 

In the next step we have included the contributions of the condensate back- 
ground self-energies in the isoscalar-scalar and isoscalar-vector channels. We 
have considered corrections possibly related to condensates of higher dimen- 
sion as well as 3-body forces, i.e. terms of order in the self-energies. The 
isoscalar condensate background self-energies turn out to be important also for 
neutron matter. However, details of the isovector interactions (beyond chiral 
one- and two-pion exchange) must still be improved. 

The relativistic nuclear point-coupling model, constrained by in-medium QCD 
sum rules and chiral perturbation theory, has also been tested in self-consistent 
calculations of ground state properties of finite nuclei. In the present work we 
have only considered spherical even-even light and medium-heavy N ^ Z nu- 
clei. Four parameters determine the nuclear matter EOS: the cut-off parame- 
ter A of the in-medium CHPT calculation of two-pion exchange diagrams, the 
two parameters of the isoscalar condensate background nucleon self-energies 
linear in the corresponding densities: G^p and Gy , and the strength param- 
eter (yfy^ of the quadratic term of the isoscalar vector self-energy. For finite 
nuclei we adjust in addition the strength parameter of the isoscalar-scalar 
derivative interaction term in the Lagrangian. This brings the total number 
of adjustable parameters to five. It has been emphasized that the empirical 
data set oi N ^ Z finite nuclei can only constrain five or six parameters in 
the general expansion of an effective Lagrangian in powers of the fields and 
their derivatives. In particular, only one parameter of the derivative terms 
can be determined by the binding energies and radii of spherical nuclei. The 
G^y turn out to be so remarkably close to leading order QCD sum rule esti- 
mates that one could have guessed their values right from the beginning. The 
remaining three parameters (the high momentum scale refiecting unresolved 
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short- distance physics, a "three body" contribution to the nucleon self-energy 
and a surface (derivative) term) behave "naturally" according to the power- 
counting hierarchy of effective field theory. Of course, questions concerning 
systematic convergence of the in-medium chiral loop expansion still remain 
and need to be further explored. 

An important result of the present analysis is that chiral (two-pion exchange) 
fluctuations are consistently included in the nuclear binding and in the sat- 
uration mechanism. These fluctuations are superimposed on the condensate 
background scalar and vector fields of equal magnitude and opposite sign, 
induced by changes of the QCD vacuum in the presence of baryonic matter. 
They drive the large spin-orbit splitting in finite nuclei. The calculated neutron 
and proton single-particle energies in ^^O and ^°Ca are in excellent agreement 
with the empirical single-nucleon levels in the vicinity of the Fermi surface, 
including the empirical energy differences between spin-orbit partner states. 

This work has demonstrated that an approach to nuclear dynamics, con- 
strained by the chiral symmetry breaking pattern and the condensate structure 
of low-energy QCD, can describe symmetric and asymmetric nuclear matter, 
as well as properties of finite nuclei, at a quantitative level comparable with 
phenomenological relativistic and non-relativistic mean-field models. This is 
a promising result, which opens perspectives of bridging the gap between ba- 
sic features of low-energy, non perturbative QCD and the rich nuclear phe- 
nomenology. 
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Table 1 

The coefficients of the expansion (40) up to order of the in-medium CHPT 
isoscalar- vector, isoscalar-scalar, isovector- vector, and isovector-scalar nucleon self- 
energies. 





V- VERTEX 


S- VERTEX 


TV- VERTEX 


TS- VERTEX 


C30 


3.96669 


3.97078 


-4.00813 


-4.12669 


C31 


-7.55255 


-7.55255 


6.54555 


6.54555 


C32 


0.762611 


0.60693 


-0.544615 


-1.73796 


C3L 


0.29288 


0.267627 


-0.370303 


-0.563877 


C40 


8.42137 


8.11713 


-2.4575 


-2.97905 


C50 


10.2384 


12.0729 


-26.8119 


-23.591 


C5L 


1.16164 


1.92167 


-7.81327 


-6.50123 



Table 2 

The coefficients of the expansion of the in-medium CHPT self energies (41) ~ (44) 
in powers of the baryon density, for the value of the cut-off parameter A = 646.3 
MeV. In this case we assumed N = Z, neglecting isovector contributions. 



Csl 


-2.805 fm^ 


Ctsl 


1.491 fm^ 


Cs2 


2.738 fm^ 


Cts2 


-1.005 fm^ 


Cs3 


1.346 fm'' 


Cts3 


-1.550 fm^ 


Cvl 


-2.718 fm^ 


Ctvl 


2.249 fm^ 


Cv2 


2.841 fm^ 


Ctv2 


-0.829 fm^ 


Cv3 


1.325 fm'' 


Ctv3 


-1.595 fm^ 
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Table 3 

Nuclear matter saturation properties: binding energy per nucleon, saturation density 
and incompressibility. The first row corresponds to the in-medium CHPT calculation 
including one- and two-pion exchange between nucleons [15]. The EOS displayed in 
the second row (PC-dd) is obtained when the resulting CHPT nucleon self-energies 
are mapped on the self-energies of the relativistic point-coupling model with density 
dependent couplings. 



MODEL 


E/A (MeV) 


Psat (fm~2) 


Kq (MeV) 


CHPT [15] 
PC-dd 


-15.26 
-14.51 


0.178 
0.175 


255 
235 



Table 4 

Nuclear matter saturation properties calculated in the relativistic point-coupling 
model constrained by in-medium QCD sum rules and chiral perturbation theory. 
In addition to the chiral one- and two-pion exchange contribution to the density 
dependence of the coupling parameters, the nuclear matter EOS shown in the first 
row includes also the isoscalar condensate background nucleon self-energies linear 
in the corresponding densities. The EOS displayed in the second row is calculated 
by including also the non- linear "3-body" contribution (proportional to p^) in the 
isoscalar vector condensate self-energy. 





E/A (MeV) 


Psat (fm"3) 


Kq (MeV) 


M*/M 


LINEAR 


-15.97 


0.148 


283 


0.753 


NON-LINEAR 


-15.76 


0.151 


332 


0.620 
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Table 5 

Binding energies per nucleon and charge radii of light and medium-heavy nuclei, 
calculated in the relativistic point-coupling model constrained by in-medium QCD 
sum rules and chiral perturbation theory, are compared with experimental values. 







E/A (MeV) 


r^^P (fm-3) 


Tc (fm 


160 


7.976 


8.019 


2.730 


2.724 


40 Ca 


8.551 


8.494 


3.485 


3.447 


42Ca 


8.617 


8.482 


3.513 


3.450 


48 Ca 


8.666 


8.396 


3.484 


3.476 


42^- 


8.260 


8.127 




3.536 


50Ti 


8.756 


8.468 




3.549 


52 Cr 


8.776 


8.494 


3.647 


3.614 


58Ni 


8.732 


8.445 


3.783 


3.748 


64Ni 


8.777 


8.451 


3.868 


3.836 


88sr 


8.733 


8.439 


4.206 


4.205 


90 Zr 


8.710 


8.443 


4.272 


4.252 
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Fig. 1. Binding energy per nucleon for symmetric nuclear matter as a function of 
baryon density. The solid curve (KFW-02) is the EOS calculated in Ref. [15] by us- 
ing in-medium CHPT for one- and two-pion exchange between nucleons. The EOS 
displayed by the dotted curve (PC-dd) is obtained when the resulting CHPT nu- 
cleon self-energies are mapped on the self-energies of the relativistic point-coupling 
model with density dependent couplings. FP-81 denotes the microscopic many-body 
nuclear matter equation of state of Friedman and Pandharipande [29], shown for 
comparison. 
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Fig. 2. Equation of state of symmetric nuclear matter calculated in the relativis- 
tic point-coupling model with and without the rearrangement contribution to the 
isoscalar vector self-energy. 
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Fig. 3. Binding energy per nucleon for symmetric nuclear matter as a function of 
the baryon density, calculated from chiral one- and two-pion exchange between nu- 
cleons (case 1, solid curve), by adding the isoscalar condensate background nucleon 
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Fig. 4. 5*2 (yo) coefficient (55) of the quadratic term of the energy per particle of 
asymmetric nuclear matter. The three curves correspond to calculations with only 
chiral one- and two-pion exchange between nucleons, with the inclusion of leading 
order isoscalar condensate background nucleon self-energies, and finally with the 
contribution of the non-linear vector condensate nucleon self-energy. Cases 2 and 
3 are explained in the caption of Fig. 3 and ZBL-99 denotes the asymmetry en- 
ergy curve obtained in the non-relativistic Brueckner-Hartree-Fock calculation of 
Ref. [36]. 
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Fig. 5. Energy per particle of neutron matter as a function of the neutron density. 
The results are obtained from chiral one- and two-pion exchange between nucleons 
(solid curve), by adding the isoscalar condensate background nucleon self-energies 
linear in the corresponding densities (case 2, dotted curve), and finally by including 
also the non-linear contribution to the isoscalar vector condensate self-energy (case 
3, dashed curve). Cases 2 and 3 are the same as in the caption of Fig. 3. This is 
shown in comparison with the microscopic many-body neutron matter equation of 
state of Friedman and Pandharipande [29] . 
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Fig. 6. Comparison of the /cj-dependence of the isoscalar vector and scalar 
self-energies resulting from a Dirac-Brueckner G-matrix calculation (solid lines: 
DB [38]), with the self-energies generated from in-medium chiral perturbation theory 
(dashed lines: CHPT + point-couplings up to 3- loop order in the energy density). 
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Fig. 7. Neutron and proton single-particle levels in ^^O and ^°Ca calculated in the 
relativistic point-coupling model. The density dependent coupling strengths include 
only the contribution from chiral one- and two-pion exchange between nucleons. 
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Fig. 8. The neutron and proton single-particle levels in ^^O calculated in the rela- 
tivistic point-coupling model, are shown in comparison with experimental levels. The 
calculation is performed by including both the contributions of chiral pion-nucleon 
exchange and of the isoscalar condensate self-energies in the density dependent cou- 
pling strengths. 
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Fig. 9. Same as in Fig. 8, but for ^Ca. 
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Fig. 10. Charge form factor (left panel) and charge density distribution (right panel) 
of ^^O. The results obtained in the relativistic point-coupling model, constrained 
by in-medium QCD sum rules and chiral perturbation theory, are compared with 
experimental data [47]. 
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Fig. 12. Neutron and proton single-particle spectra of ^"Ni, calculated with the 
standard relativistic mean-field model NL3 effective interaction [48], and with the 
relativistic point-coupling model constrained by in-medium QCD sum rules and 
chiral perturbation theory (FKVW). 
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